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Restructuring of electronic spectrum in a buckled silicene monolayer under some applied voltage 
between its two sublattices and in presence of certain impurity atoms is considered. A special 
attention is given to formation of localized impurity levels within the band gap and the to their 
collectivization at finite impurity concentration. It is shown that a qualitative restructuring of 
quasiparticle spectrum within the initial band gap and then specific metal-insulator phase transitions 
are possible for such disordered system and can be effectively controlled by variation of the electric 
held bias at given impurity perturbation potential and concentration. Since these effects are expected 
at low impurity concentrations but at not too low temperatures, they can be promising for practical 
applications in nanoelectronics devices. 
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I. INTRODUCTION 

After revolutionary breakthrough of graphene, intro¬ 
ducing unusual relativistic effects into solid state physics 
[ij, the family of relative materials is continuously grow¬ 
ing. In particular, a possibility of obtaining new semi¬ 
conducting materials where the bandgap can be tuned 
by external electric bias is extensively studied. This was 
first demonstrated for the graphene bilayer, called bi¬ 
graphene 0. Unlike the basic monolayer graphene, here 
a non-equivalence of two sublayers takes place under elec¬ 
tric field applied normally to them. Opening of a tun¬ 
able semiconducting gap and perspectives of its practi¬ 
cal use in tunable transistors is now broadly discussed. 
The same possibility was already indicated in a variety 
of similar systems, including even single layered, as, for 
instance, silicene, the Si-based analog of graphene 3- Its 
important structural difference from graphene consists in 
a much more pronounced buckling of its 2D hexagonal 
lattice, thus leading to a non-equivalence of two sublat¬ 
tices in the same layer and to opening of a bandgap un¬ 
der normal-to-plane electric biasj3-@, Subsequently, 
fabrication of practical field-effect transistors based on a 
silicene sheet is extensively sought for 01|- 

When comparing these 2D systems with common semi¬ 
conductors, an important question arises on their behav¬ 
ior under doping by impurity atoms. As well known, 
such doping in common semiconductors produces local¬ 
ized in-gap energy levels near the edges of conduction 
band (called donor levels) or valence band (acceptor lev¬ 
els) [i^. For mostly used dopants (as Si neighbors from 
the periodic table) these levels are very shallow (of some 
tens meV depth compared to some eV bandwidths) so 
that charge carriers can be thermally excited from them 
to the nearby band (conduction or valence) and thus con¬ 
tribute into conduction of respective kind (electron or 
hole). The resulting conductivity turns sensitive to ex¬ 
ternal bias realized in specific devices, defining their ef¬ 
fectiveness Ha. Typical dopant concentration n is quite 
low, in order to assure the mean distance f = 


between dopants to surpass the long effective radius of 
localized state rioc ^ a (the lattice parameter), it should 
not exceed hq = ^ 10“^^ cm“^ (for 3D systems). 

Then it is known that for n ^ uq the doped system 
is brought to metallization [l^ . due to growing interac¬ 
tion between dopants and subsequent broadening of the 
dopant level. Hence the Fermi level, initially fixed at 
the dopant level, becomes displaced to the band inte¬ 
rior. This phenomenon is generally considered adverse 
for electronics purposes, since it drastically reduces the 
bias sensitivity. Otherwise, an alternative type of impu¬ 
rities (as transition and rare-earth elements), producing 
the so called deep levels in the semiconducting gap [l3l |. 
are not effective for thermalization of carriers and instead 
can act as traps for them. 

The above limitations however can be effectively over¬ 
come under the possibility for tuning the fundamental 
bandgap and also other relevant spectrum characteris¬ 
tics as the Fermi level and the Mott’s mobility edges M- 
This opens a formerly unexplored perspective of bring¬ 
ing the bias sensitivity of conduction to a much broader 
scale then in common semiconductors, spanning it from 
metallic to insulating regimes through the Mott’s metal- 
insulator transition (MIT). The purpose of the present 
study is to illustrate such an expectation on the partic¬ 
ular example of biased and doped silicene, considering 
both situations of shallow and deep dopant levels with 
their specific regimes. Though a detailed treatment of 
these issues, using the realistic impurity potentials, their 
screening by relativistic electrons, etc., can present cer¬ 
tain technical problems, it can be much facilitated with 
use of simplified models, traditionally applied for impu¬ 
rities both in common semiconductors and in grap hene- 
related materials. Such are the Lifshitz model [I 3 , bet¬ 
ter suited for shallow dopants, and the Anderson hybrid 
model 0, more adequate for deep dopants. Below we 
analyze the electronic spectra of biased and doped sil¬ 
icene within these two models and indicate possible tun¬ 
ing regimes to reach desirable effects. 
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FIG. 1: Crystalline structure of a buckled silicene plane where 
silicon atoms are up- or down-shifted hy 12 from the initial 
hexagonal plane (red lines). The dashed line delimits a unit 
cell at the position n with two non-equivalent sites, of A-type 
(clear) and B-type (dark), blue arrows indicate three vectors 
8 between nearest neighbor Si atoms. 


II. FORMULATION OF THE PROBLEM 

For the silicene hexagonal lattice with two non¬ 
equivalent sites in unit cell (Fig. [T]), we write down the 
tight-binding Hamiltonian as: 

-ffo = (1) 

k 

Here 2-spinors ipl. = 6]^^ are made of 2D 

Fourier transforms Ok = e*’^ "an and 6k = 

e*’^'"6n of local Fermi operators at A- and 
B-type sites in the n-th unit cell {N the total num¬ 
ber of cells). The 2x2 Pauli matrix expansion hk = 
(T 3 H /2 -I- i 7 +tk + includes the on-site energy shifts 
±H/2, due to the effect of buckling and external electric 
field, referred to as bias in what follows, and the com¬ 
plex factors tk = ‘ikie to the hopping ampli¬ 

tude t ^ 1.1 eV between nearest neighbor sites separated 
by the vectors {S,dz) with their xy-plane components 
1^1 = a ~ 0.22 nm and normal buckling components 
dz ^ 0.2a [13) [3 • For sake of simplicity, relatively weak 
spin-orbit interactions and spin degrees of freedom are 
omitted in Eq. [T] The relevant low-energy physics is gen¬ 
erated near the nodal points ±K = (±47r/3-\/3a, 0) in the 
Brillouin zone so that for q = k — K with ag <C 1 we have 
tk ~ where the Fermi velocity vp = 3at/{2h) 

and (pq = arctangy/g^,. 

We study electronic states in this system using the 
Green function (GF) matrix G(k, k') = (('0k|'*/'k'))that 
satisfies the equation of motion: 

eG(k,k') = (2) 

for the full Hamiltonian H. Generally, the GF matrix 
defines the system energy spectrum by the roots of the 
general dispersion equation: RedetG = 0, that is by the 


poles of its diagonal elements in the basis of exact eigen¬ 
states, and the total density of states (DOS) is defined 
as: 

P(e) = ^;^TrImG, (3) 

in any basis. Using of Eq. [5]with the unperturbed Hamil¬ 
tonian, H = Hq, Eq. [H leads to a momentum-diagonal 
form for the non-perturbed GE matrix G^'’^(k, k') = 

(5k,k'Gk°^ where: 


MO) ^ s + hj^ 

^2 _ ^ 2/4 _ ^2 > 


(4) 


and ^k = l^kl is an effective momentum variable for 
the relevant low-energy range, jej <C A (where A = 

hvp y^Kja = t\/Tr\/3 is the bandwidth, that is an energy 
cut-off). In this approximation, the electronic spectrum 
consists of two symmetric bands ±eo (^) = ±-\/U^/4 -|- 
with a gap of V between them. The related non- 
perturbed DOS is linear beyond this gap: 

po (e) = (UV4 - 6^) 9 (A^ + UV4 - e ^), (5) 


and normalized: Po (s) de = 1. This situation in 

silicene with buckled lattice is referred to as realization 
of a tunable gap between the lower (valence) band and 
the higher (conduction) band, in contrast to a fixed gap 
in common semiconductors (like crystalline Si). 

Now we pass to the specifics of doping this system by 
impurity atoms at random sites within its crystalline lat¬ 
tice that produce its perturbation through the above re¬ 
ferred Lifshitz and Anderson models. 


III. LIFSHITZ MODEL 

We begin with the simpler Lifshitz model [3 charac¬ 
terized by a single perturbation parameter, the on-site 
energy shift U on all impurity atoms located at pi sites 
in the 1st sublattice and p 2 sites in the 2nd sublattice 
with relative concentrations ci _2 (both expectedly small, 
Cj <C 1). Such substitutional type better corresponds to 
impurities such as Si neighbors in the periodic table. The 
corresponding Lifshitz perturbation reads: 

(6) 

k,k' j,pj 

with the scattering matrices Uj = Upj and the sublattice 
projectors pip = (1±(Tz)/2. In this model with the 
full Hamiltonian H = Hq + the explicit equation of 
motion for the momentum-diagonal GF matrix: 

G(k) = G(o)(k) 

+ ^ E e*('^'-‘^)p^4‘’^G,G(k',k), (7) 
















3 


leads to the standard solution: 


Gk 


~ ^l,kPl — S2,kP2 


( 8 ) 


Here the partial self-energy functions are presented by 
their respective group expansions (GE’s) 


'S'i.k — CjTj (1 -|- CjBj ]^ + ■■■)> (9) 


where the partial T-matrix Tj = U/ (1 — UQj) with lo¬ 
cal GF’s gj = (Gk) describes the effects of 

c P 

multiple scatterings on single impurity center. It can be 
generally shown that the GE series, Eq. [HI is converg¬ 
ing within the energy range of band-like states and can 
be well approximated there by its first T-matrix term 
while the rest of terms are important for the check of 
convergence. The first of non-trivial GE terms is due to 
impurity pairs: 


^i,k 




—zk-] 


“1“ —n 


n/0 


1 


( 10 ) 


and includes the functions Pj,ii = 

^k'/k " (Gk'^ of inter-impurity inter¬ 

action. It should be also noted that all the products of 
these functions in the expansion of Eq. [TU] are presented 
by multiple sums in non-coinciding momenta ( 2 ^, [ 2 l|. 
The omitted terms in the brackets of Eq. |9] correspond to 
clusters of three and more impurities, they are expressed 
through respective combinations of these functions. 

Generally, once the initial translation symmetry of 
crystalline lattice is broken by the presence of impuri¬ 
ties, the quasi-momentum is no longer an exact quan¬ 
tum number and also the system spectrum is not lim¬ 
ited to the initial bands, since already a single impurity 
can produce localized levels beyond the bands. However, 
for not too strong a disorder, this spectrum maintains 
continuous ranges of band-like states (both the modified 
initial bands and possibly some new, impurity, bands 
arisen near localized levels), intercalated by the ranges 
of truly localized states (either on single impurities or 
on their clusters). Within band-like ranges, the approx¬ 
imate dispersion laws of corresponding subbands Sj (^) 
(in our case, j = 1,2, imp for two initial and impurity 
bands respectively) are given by the formal roots of the 
above mentioned dispersion equation in the k-basis: 


RedetGk = 0. (11) 


ending validity of the IRM criterion. This also qualita¬ 
tively agrees with a similar change in the convergence 
criterion for GE: c|i?yk| ^ 1- Within the bandgap, the 
broadening Tj mainly results from ImHj k in Eq. so 
analysis of this range needs calculation of the functions 
gj and „. A reasonable approximation for them fol¬ 
lows from substitution of Gk by g|^°^ in corresponding 
sums resulting in: 


91,2 


e±Vl2 


In 



(13) 


and: 




(.e±V/2)T, 

IP 


Kq {n/rPj 


(14) 


(see details in Appendix). Here the characteristic length 
Tg = hv-p / !A — and the McDonald function 

ATp {x) has asymptotics 


ATo {x) 


In ( 2 /a:) — 7 , a; <C 1 , 
a/2/ ( 7 ra:)e a: > 1, 


with the Euler’s constant 7 Ri 0.5772. 

The logarithmic divergence of gj, Eq. [T31 near one 
of the gap edges allows a localized level eioc to appear 
there under a proper impurity perturbation. Thus, if 
one chooses for definiteness U < 0 (and supposedly |17| < 
A), this level is due to the pole of Ti (by impurities in 
1st sublattice) near the upper edge E/2, their separation 
being well approximated as: 

r/2-£,„«.Ae-V/(|K|v-)sc„/. (15) 

Hence the localized level is exponentially shallow for all 
practically achievable bias values (always E <C A), which 
can justify such modeling of real shallow levels. At this 
choice, another term T 2 (by impurities in 2 nd sublattice) 
has no poles and is less relevant. 

Further analytical study of modified spectrum uses 
some approximated energy dependencies of the relevant 
Ti-matrix. Thus, in a close enough vicinity to the local¬ 
ized level, |e — eiod -C E/2 — Sioc, its denominator can be 
linearized: 


A^ (E/2 - 
V{e- £ioc) 


(16) 


However their validity is restricted by the known loffe- 
Regel-Mott (IRM) criterion [T3 . that for this case 
takes the form: 

^de,(C)/d^»r,(0 (12) 

with the damping term Tj (^) = ImA^ k for ^ = ^k- The 
separation points between extended and localized ranges, 
called Mott’s mobility edges, are estimated from the con¬ 
dition that the symbol “<C” in Eq. [T^is changed for 


while in a wider area, 0 < E/2 — e <C E, the logarithmic 
approximation applies: 


Ti 


A2 _^vl2-eioc 

— In -^-. 

E E/2 - e 


(17) 


At last, when considering the energy scales over the whole 
band gap, |e| > E, the complete formula, Eq. [T31 should 
be used in the T-matrices. 
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FIG. 2: Spectrum restructuring in silicene with Lifshitz impurities at the choice of their parameters U = —A and c = 0.01 in 
function of the bias V. Insets show the particular dispersion laws at: a) weak {V = 0.004A), b) medium {V — O.llA), and 
c) strong {V — 0.23A) bias. Blue areas present filled band-like states and dark yellow do filled localized states, red dashed 
lines mark the concentrational broadening of impurity level. Double dashed lines indicate the critical bias levels Vdec for band 
decomposition (orange), Vm for Mott (blue), and Va Anderson (purple) phase transitions. 


Of course, for t7 > 0, symmetric formulas with respect 
to Eqs. [IHl mi [TT] apply for the pole of T 2 near the lower 
gap edge and for its related vicinities while Ti becomes 
irrelevant. 

Let us focus now on the most restructured region of 
spectrum, including the impurity band SimpiO s-nd its 
closest neighbor areas of £ 1 , 2(0 bands. In this course, it 
is convenient to consider this restructuring with growing 
bias V at fixed impurity parameters U and c. The T- 
matrix approximationZ'j k Ri CjTj is sufficient at the first 
step, as far as the quasiparticle lifetime and respective 
IRM limits for band-like states are not considered. 

At lowest bias V <C c|17|, the numerical solution of the 
dispersion equation, Eq. [TTl with use of Eq. mi shows 
the spectrum restructuring to be very close to its sim¬ 
ple shift by cU, that is the impurity effect is reduced to 
that of effective medium potential. Formally, this solu¬ 
tion includes the impurity band £imp(0 of o|C/| width 
and separately the modified upper subband £i(0' How¬ 
ever, it is seen from Fig. [^Kinset a) that their compo¬ 


sition is closely matched near the Sioc level and prac¬ 
tically coincides with a single shifted law: Scomp (C) = 
i/E^/4 -I- -I- cU, attaining its lower edge at ^ = 0: 

£g = £comp{0) = H/2 + cU. As seen from Fig. [2jinset 
b), such composite £1 -|- £imp-band structure persists for 
low enough bias such that c <C Cq, and, using Eq. 1151 
this relates to: E <C lA = A^/ (|C71nc|)(the latter value 
to be explained below). Noting from Eq. [141 at £ = £ioc 
that the localization radius is rioc = fivp / (Ay^), this 
refers to rioc r = hyp/ (AyAc), a 2D analogy to the 
metallization condition discussed in Introduction. Thus, 
the characteristic concentration cq (in fact, an analog to 
riQ in Introduction) is bias tuned, and a tuned Anderson 
transition on the £imp-band (mer ger of its mobility edges 
and vanishing of band-like states ) is expected at the 
bias value Va- 

At growing bias to V ^ c|C/|, the composite band 
structure can be described with use of Eq. [T71 Then 
the renormalized dispersion law £ = £imp (C) follows from 
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Eq. [TT]in an implicit form as: 



This equation permits analytic solutions near the edges 
of Eimp-band. Thus, the lower edge corresponds to ^ ^ 0: 


.( 0 ) 


V cA2 
J~ VW (c/co) ’ 


(19) 


where W (z) is the Lambert W-function [ 2 J. Its asymp¬ 
totics, IT (z » 1) ~ In (z/ Inz), used in Eq. [T^ provide a 
simpler function: 


Using it in the IRM criterion, Eq. [T^ we estimate the 
location of the mobility edge £c near the bottom of £imp' 
£c — £g ct/^U/A^. It is negligible beside the gap of V 
between the £ 2 - and ^imp-bands and the width Aimp of 
£imp-band. However, the broadening, Eq. [531 at e ~ £ioc 
much exceeds the formal gap ~ Vl^ — £ioc (here exponen¬ 
tially small) between £imp- and ei-bands. This permits 
to consider such gap and the very level £ioc non-existing 
and justifies the concept of a composite band in the weak 
bias regime (see insets a,b in Fig. [5]). 

The overall electronic state of the doped system is de¬ 
termined by the location of its Fermi level ep with respect 
to the mobility edges. Supposing each impurity atom to 
supply one carrier to the system and its undoped state 
to possess eF=0, its position at finite c is found from the 
equation: 


V Tj I In c U 
2 . In (cco) Va 


( 20 ) 



(25) 


replacing the above linear £g dependence at U <C c|?7| in 
the wider range of U <C Ua- At yet higher bias, up to 
V ^Va, the full Eq. [T51 applies. 

Expanding Eq. [TS] at <C U^/4 — e^, the long-wave 
dispersion law is obtained: 


-imp 


(0 


Eg -I- I 1 


y 

4cA2 j 




indicating the ecomp-behavior for V <C \/cA. The further 
growth of £impiO finally reaches the short-wave asymp¬ 
totics (at » cA^): 


£^n^p(0-£^oc-cA" ^^^^/'°^ (22) 

This defines its formal upper edge £f = £imp{A) •££ £ioc — 
c{Vj2 — £ioc) [13, and then the total width of impurity 
band: 


A,mp-£f y (^(c/co) (23) 

Growing with V from the initial value of c|17|, A^^p 
by Eq. |23| would reach a maximum at some U* = 
E(c)A^/|{7| with the factor F (c) varying from m 0.12 
to Ri 0.22 in the range of 10“"^ < c < 0.1. However, such 
14 is already close to the critical value Va and hence to 
the impurity band collapse, making this maximum mean¬ 
ingless. 

The next step is to determine the lifetimes of the ob¬ 
tained quasiparticle states, in order to establish the IRM 
limits for their existence. For the quasiparticle with en¬ 
ergy £ = £comp (C), <3an consider its effective broaden¬ 
ing: 


(including the spin factor 2). For the weak bias regime 
(or, in other words, for cq -C c), one can use here the 
composite band DOS: p (e) ~ po (£ + £g), and obtain the 
bias dependent Fermi level within this band as: 


£f (U) r: y^cA^ -I- £2. (26) 


At low enough bias, V <C v^A, it lies as high within the 
Ecomp-band as £f R v/cA (Fig. j^) and defines a metallic 
behavior of the system. This can be just compared to 
metallization of common semiconductors at high enough 
doping (c » Co). 

With growing bias, £f (V) gets closer to the band’s bot¬ 
tom £g but its expected crossing of a mobility edge and 
the system transition into insulating state can be only 
reached after the decomposition of £imp and £1 bands by 
means of an emerging localized range around the impu¬ 
rity level £/oc- From comparison of Eqs. [T5land[24l this is 
estimated to take place at c ~ cq In^ (I/cq). The related 
bias value is high enough: 


VtZec 


A2 

2|UW_i(-yB/2)|’ 


(27) 


including the lower branch W_i of the multivalued Lam¬ 
bert function with asymptotics lU_i (z) R ln(z/ln|z|) 
for—l/e<z<0 ( 23 . However, this Vdec is yet well be¬ 
low the £ij„p-band collapse value Va (see inset c in Figj^)- 
Then, taking into account the T-matrix contribution to 
DOS for the lower decomposed subband : 


Pi (e) « 


£ 

A2 


CCqA^ 

2U" (£ - £loc) ’ 


(28) 


and using it in Eq. 1251 we find the condition that £f 
reaches the top of £imp' 


comp 


(£) = Im 


cU 


1 - Ugi (£ - cU) 

UV {V/2 + £-cUf 
Ay £ — cU 


(24) 


r^f 

C = 2 Pi{£)d£ I 

Jsg 

CCqA^ 


+ 


■In 


2 U 2 V CO IT (c/co) 


IT (c/co) 
1 1 
c 


- Co 


TTC 

Y 


(29) 


















Implicitly, Eq. [551 defines the bias value Vm that can be 
associated with the tuned MIT, provided this value be 
above Vdec so that the top of Simp already pertain to the 
localized range. The MIT bias value is estimated from 
Eq. [551 as: 


Vm 


7A2 

4|t/| In (xm/c)’ 


(30) 


where the factor in the logarithm depends on the pertur¬ 
bation parameter as: Zm ~ (3.35C//A)^, by a reasonable 
empirical fit. Then, the numerical comparison between 
Eqs. [571 and 1551 shows that Vdec in fact precedes Vm for 
all realistic U < A. But the sequence of tuned MIT and 
Anderson transitions can be changed depending on the 
impurity parameters. So, the Va —^ Vm sequence for 
their above choice (as in Fig. HD passes to Vm —t Va 
for c = 0.01 and \U\ > 0.71A. With the bias V exceed¬ 
ing Va, the impurity band does not make sense already 
but there exists a well defined localized level, Eq. [151 
whose width at c <C cq estimated from the GE pair term 
(see details in Appendix) becomes exponentially small: 

fooc - Co (AVe) 

The above discussed tuned restructuring of spectrum 
can be summarized as follows. At low bias, V <C Vm, 
the system is metallized by the impurity doping, with 
the Fermi level lying deep within the composite ei+eimp- 
conduction band. The composite band gets split into ei- 
and e^mp-bands, separated by the range of localized states 
around the impurity level Sioc, at bias reaching Vdec- At 
further growing bias to Vm, the Fermi level meets the 
mobility edge above sioc to produce MIT. After the An¬ 
derson transition occured at E = Va, a single £i conduc¬ 
tion (unoccupied) band is left in the spectrum, the Fermi 
level staying fixed near eioc- All the critical bias val¬ 
ues, Vdec, Vm, and Va, can be reduced by choosing lower 
impurity concentration, but this reduction is as slow as 
^ 1 / In (1 /c) and simultaneously the thermal stability 
level for MIT is reduced as k^Tmax cln (1/c) A^/|[/|. 

Now let us consider the alternative scenario, or the 
Anderson hybrid model. 


IV. ANDERSON HYBRID MODEL 


In the Anderson model, there are two perturbation pa¬ 
rameters: the on-site energy eg for an electron at an im¬ 
purity atom and its modified hopping amplitude Tyf (sup¬ 
posedly ?7 ^ 1) to the nearest neighbor host sites. Such 
type of impurity perturbation better corresponds to tran¬ 
sition metal or rare-earth atoms, known to produce deep 
levels in common Si. In silicene, these atoms predomi¬ 
nantly occupy interstitial positions, linked to both host 
sublattices (see Fig. o. In the corresponding perturba¬ 



FIG. 3: Anderson impurity at an interstitial position in sil¬ 
icene lattice. 


tion Hamiltonian: 


Ha 


'y [sockpCKp 

p 


1 


^ + h.c.) 


(31) 


this linkage is presented through the 2-spinor = 
The other difference of this model consists in 
the presence of independent Fermi operators ctp for an 
electron on impurity site p, besides the above considered 
host operators in the '0k spinor. Subsequently, it gen¬ 
erates more involved GF structures (with respect to the 
sublattice indices): besides the “host” ((((k, k') matrices 
by Eq. El here also “impurity” gp.p' = ((aplctp/)) scalars 

and “mixed” /ik.p = ((V^klaj,)) and = ((Q;p|0^)) 

spinors appear. 

Under the full Hamiltonian H = Hq-\-Ha, the equation 
of motion for the “host” GF matrix gets modified from 
Eq. [7] to: 


G(k,k') =4.k' G(°>(k) 

^ P 

and the respective equation for the “mixed” spinor in 
its right hand side reads as: 

hU' (^ - ^o) = ^ E e*‘^" PTt„G(k", k'). (33) 

VA 

Then the specific solutions for all the above mentioned 
GF types follow from consequent iterations of Eqs. [35[ 
andj^ 

The general strategy consists in separating, after each 
iteration step, all the GF’s in the right hand side that 
were already present in the previous steps and doing 
next iteration for the resting ones, in order to com¬ 
pose and then fully solve an equation for each GF. The 
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most important between them are the diagonal G(k) and 
p that enter the total DOS 

in the generalization of Eq. |3l Thus, a full solution for 
G(k) follows from Eq. [35] (at k = k'), after separating 
the same G(k) in the right hand side of Eq. [33] for Zip 
and applying again Eq. |32| to all G(k',k) with k' 7 ^ k 
there. At this next iteration, a similar separation of Zip 
is also done, giving rise to a respective full solution for 
Zip and so on. The result, formally analogous to Eq. jS] 


G(k) = 



(35) 


includes the self-energy matrix in the GE form, similar to 
Eq. m .Ek = cTk (1 -I- ci3k -I- ...). But here the T-matrix 
term: 


,2f2 


cTk = 


V i 


1 + 2772 IV ^ e 

p 


1 ^ l-t-d+e^**^ 


tT_e 


—2i(p 


^imp ^^imp 


(36) 


has an important difference from the Lifshitz model 
form, Eq. |9l in its momentum dependence, both on 
the radial variable ^ = ^k and on the angular argument 
ip = (fk-K- It also includes the single impurity level 
£imp = £ 0 / (1 + 2 ? 7 ^) (reduced by its coupling to the host) 
with its imaginary part Timp = Im and the 

“impurity” scalar self-energy: 


Ep — ^ Ap_p/ [E'p/„p 

p't^p 

+ ^ ^ Ap/_p// (E'p//_p -|-...) 

p"#p,p' 


(37) 


Here the scalar functions 
1 

N 


F, 


p-p 


Y.eMp-p')rlG^r^, 


A, 


F 


p-p' 


£ ^imp ‘^^imp Ep/ 


describe the effects of indirect interactions between impu¬ 
rity centers. The latter Ap_p' functions also define the 
GE terms of the “host” self-energy Ak, along the same 
formal structure as in Eq. da while the GE structure for 
the scalar Sp in Eq. |37|is notably different. In this way, 
the solution for the diagonal “impurity” GE follows as: 

~ (l + 27?2)(e-e,^p-ir,™p-Ep)’ 

and can be then used in Eq. [M] A specific feature of Ep 
is the random statistical distribution of its values due to 


random p' positions around given p, with the standard 
deviation a — y^E^ — Ep^. Since a finite range of inter¬ 
impurity interactions, this deviation does not vanish in 
the thermodynamical limit ZV —>■ 00 , unlike that for Ak 
(known as the self-averaging property HI). As to the 
mean self-energy Ep, it can be simply included into the 
definition of impurity level Si^p, so understood in what 
follows. 

In this model, we present the basic secular determinant 

as det (^Gk^ = det ^Gk°^^ -I-Ek, with the scalar self¬ 

energy: 


Ek = 


(£ + ?cost/?) 


N 


1 + 2772 

E- 


£ ^imp ZT/^p Ep 


(39) 


Then the general Eq. dl (written at the T-matrix level 
in neglect of T/mp and Ep) reads: 


= + (40, 

£ ^imp 

where the reduced impurity concentration c = 
erf'/ (1 + “Irf^ measures the quasiparticle weight trans¬ 
fer from impurity to band states. Eq. 001 defines the dis¬ 
persion laws for “host” quasiparticles that are quite close 
to the non-perturbed eo {^) except for the Simp vicinity 
of ^ csimp width where the splitting of two subbands 
is mainly developed. Within that splitting range, both 
subbands strongly deviate from Eq (^) and display a sen¬ 
sible in-plane anisotropy: £± (^, </?) (Fig. 0]), unlike the 
isotropically split subbands in Fig. 01 Physically, this 
anisotropy reflects the breakdown of local inversion sym¬ 
metry for an impurity interstice at applied field bias. 

Other difference from the Lifshitz model is in the possi¬ 
bility that, at varying bias E, the impurity level eimp can 
be crossed by the band edge E/2. Lastly, the presence 
of factor in the T-matrix, Eq. |36l leads to vanishing 
damping for long-wave quasiparticles, so that the mobil¬ 
ity edges should correspond to shorter wavelengths (to be 
defined below). This implies that Bloch-like states with 
such short wavelengths would not make sense at energies 
close enough to the impurity level. 

As a result, the composition of electronic spectrum in 
the Anderson model is more complicated than in the Lif¬ 
shitz model. Here we have generally up to three subbands 
of the states by electrons on host sites: the two split 
£± (Cl +) subbands and the almost non-perturbed —Eq (C) 
(valence) subband, together with a subband of the states 
on impurity interstices (Eq. 1381) . The important sepa¬ 
ration between band-like and localized states can be es¬ 
tablished from the following principle. A given energy 
E pertains to the range of band-like states by virtue of 
those solutions of Eq. |33| that satisfy the IRM criterion, 
and if no such solutions exist, this energy pertains to the 
range of localized states by virtue of the related solutions 















Now, to evaluate the standard deviation u, we restrict 
Sp, Eq. I37L to its first term and approximate the in¬ 
teraction function by using the non-perturbed GF (see 
details in Appendix): 




X 


N ^ 


k 


/ n 

\ . ncos9 

eKo - 

1 

> 

[ Vu 

/ 



(43) 


Since within the relevant energy range for this case, 
E/2 < e, the argument of McDonald functions turns to be 
imaginary, they can be exp ressed through the 1st and 2nd 
kind Bessel functions (^ : Kq (ix) = ^ [Eb (x) + iJo (x)] 
and Ki {ix) = f [— {x) -I- iYi (x)]. Notice the presence 
of p-wave anisotropy by cos 9 = Uxln in Eq. I43L similar 
to that of the dispersion law, Eq. IJHl Next we obtain 
(see Appendix): 


FIG. 4: Silicene dispersion laws near their splitting by 
the impurity level Simp (red dashed line) show an in-plane 
anisotropy. Inset resolves them for particular directions: 1 
for e+(^,7r), 2 for e+(C,0), 3 for e_ (^,0), 4 for -(5,7r/2), 
5 for £_(^, tt). The impurity parameters are chosen as: 
Simp = O.IA, c = 0.005 and the bias V = OAeimp- The arrows 
indicate the splitting range near the impurity level Simp- 


^ cY,iAnF-n)^ 
n^O 

{oo^imp) ^ In ^ 
(S — Simp)^ Ccr 


(44) 


of Eq. [351 All these states contribute to the total DOS 
with their particular weights defined by the residues of 
corresponding poles of diagonal GF’s in Eq. [M] 

Following this principle, one can apply the IRM crite¬ 
rion, Eq. m for almost isotropic band-like states beyond 
the splitting range, |e — Simpl cSimp, but will need its 
more complete form: 


k- Vk£± {i,^) » r± (^,p), 


(41) 


within this range. Expecting the most important mobil¬ 
ity threshold to be located below Simp, we focus on the 
relevant damping term r_ {^, ip) = Im Ek|e — eimp|/|e^ — 
E^/4|. There are several additive contributions to this 
term, due to Timp and Ep in the denominator of Tk, 
Eq. [331 and to the GE terms such as i?k- A more de¬ 
tailed analysis shows that the dominating contribution 
comes from Ep, expressing the decay rates of band quasi¬ 
particles into the randomly distributed localized levels 
£imp+Sp. This contribution is already self-averaging and 
its average involves the probability distribution function 
P(Ep). Since Ep by Eq. [33 is a sum of a great num¬ 
ber of independent random values like Ap_p/Pp/_p, this 
probability distribution should take a normal form: 


Pi^p) = 


1 




-exp 


that readily implies: 


(Ep - Ep) 


2 ct 2 


Im Xlk — (s' ~h ^ cos Im (s ^imp H” ^p) 




(£ ^imp) 

2^ 


• ( 42 ) 


presenting an energy-dependent cr (e) that grows at ap¬ 
proach to the single impurity level Simp, as it can be 
expected for such resonance interactions. In Eqs. |331 
l44l Co = \sjmp ~ E^/4|/A^ generalizes the definition in 
Eq. [T3|for the locations of band edge E/2 either above 
and below e^mp, and the meaning of Ccr = ^’Impl^ i® 
explained in what follows. Using the result by Eq. 031 
in Eq. I42L we conclude that the relevant damping term 
r_ (^, p) very steeply shoots up from exponentially low 
to values as high as ^ [cIcq) (e -f ^ cosv?) Eimp when 
reaching the condition |e — eimp\ cr (&) or: 


^ ^imp 



1/4 


^O^imp ■ 


(45) 


So Eq. 03] gives just an estimate for the distance from 
£imp to the mobility edge We notice that, for c close 
to Ccr and cq, this distance exceeds the width cSimp of 
the splitting range that justifies the above usage of the 
non-perturbed spectrum in Eq. 1431 Also, the estimate 
by Eq. 05] defines the above referred minimum admit¬ 
ted wavelength of band quasiparticles with energies near 

^imp- ^min ^ o/v^. 

The principal practical issue of tuned metal-insulator 
transition is resolved by comparing the mobility edge Sc 
and the Fermi level ep whose initial position at zero bias 
is below Simp if the impurity concentration is smaller of 
the above defined critical value: c < Ccr- 
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FIG. 5: Dispersion curves for silicene with Anderson impurities as in Fig. |4]at growing bias V that moves the bottom of the 
conduction band towards the impurity level Simp- MIT occurs at Vm ~ 0.84eimp (in between a and b panels) and the Anderson 
transition at Va ~ 2eimp (in between c and d panels). 


Expecting the crossing of ep with Sc at growing bias to 
occur outside the splitting range, we can safely estimate 
ep with use of the non-perturbed DOS, Eq. [SJ in Eq. [5S] 
to result in: 

ep (E) v'cA2 + v^/4, (46) 

instead of Eq. [521 Then the MIT bias Vm, when the 
Fermi level crosses the mobility edge, is found for the 
Anderson model as: 


Vm - 2e,^pJl -—(l + In^/^ —) , (47) 

y ^CT \ ^cr / 

that i s slightly below the critical value Vcr = 
2eimp\/^ — cjCcT, when the Fermi level reaches Eimp (as 
shown in Fig. 0. With further growing bias, V > Vcr, 
the Fermi level stays fixed near £imp while the Anderson 
transition for the (^, (p) subband takes place when the 
mobility edge is attained by the band edge E/2. This 
corresponds to Va ~ 2£imp- 

Finally, at E > Va, the impurity level £imp stays below 
the bottom of the almost unperturbed main band (as in 
Fig. mi), then Ep and so the broadening of Simp becomes 
exponentially small ^ coAe“'^“/'^ (at cq ^ c) by virtue 
of similar decay of in Eq. |33| and there are almost 
unperturbed subbands ±£o (C) ia fhe spectrum (alike the 
above case of the Lifshitz model). 

V. DISCUSSION AND CONCLUSIONS 

The above presented considerations of electronic spec¬ 
trum in biased and doped silicene within the frameworks 
of two models for impurity perturbation show a variety 
of restructuring processes in this spectrum with different 


dynamics for its particular subbands, derived from both 
the initial host subbands and from the impurity levels. 
The main difference between the two models is in the lo¬ 
cation of impurity energy level. In the Lifshitz model, it 
is Sloe, Eq. [m which closely follows the biased edge of 
one of the main subbands while in the Anderson model it 
is the bias independent £imp, Eq. 1331 that can be crossed 
by the biased main band edge. This determines the dif¬ 
ference in spectrum transformations with growing bias 
E. 

In the Lifshitz model, the initial metallic state at 
E <C Em corresponds to the Fermi level £p well above 
£/oc and the tuned MIT is realized through its dropping 
down to the mobility edge that emerges near eioc after 
the impurity subband gets decoupled from its neighbor 
main subband at E > Vdec- This process develops rather 
slowly with growing bias and requires the stronger critical 
level Em the higher impurity concentration c is present in 
the sample. In contrast, the Anderson model provides a 
possibility for initial £f to be positioned below Simp and 
to reach the mobility edge with growing bias, the sooner 
the higher c is chosen. These scenarios can be suitably 
presented in the form of phase diagrams in terms of the 
relevant variables “impurity concentration-electric bias” 
(Fig. 0. 

These diagrams for the two models clearly display the 
above mentioned difference in electronic phase dynam¬ 
ics. In the Lifshitz model, the critical value as a func¬ 
tion of impurity concentration, Em (c), grows (logarith¬ 
mically slowly) (Fig. 0i), while in the Anderson model 
this function is rapidly decreasing from its initial value 
Em (0) = 2eimp (Fig. 03). Also, there is a notable dif¬ 
ference in behavior of other phase boundary in this case 
which defines the Anderson transition in the collapsing 
impurity band. In the Lifshitz model, Va (c) grows in a 
similar way to Em (c) and, depending on the perturba- 
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promising for such purpose. Comparing the present sys¬ 
tem to the other known material with tunable gap, the 
bigraphene, where similar doping effects were recently 
considered [2^, an advantage of the siliceiie host is seen 
in the simpler structure of its electronic spectrum. Ex¬ 
perimental checks on the proposed regimes of doping and 
tunable phase transitions could better determine the field 
for future studies and probably open some new possibil¬ 
ities in this direction. 


Acknowledgements 


FIG. 6: Phase diagrams of electronic states in biased doped 
silicene in the variables “impurity concentration-electric bias” 
for two models of impurity perturbation. The blue areas cor¬ 
respond to metallic phases, separated by the Mott MIT lines 
Vm (dark blue) from insulating phases (white areas), Ander¬ 
son transitions (collapse of impurity band) are shown by the 
red lines, dashed lines indicate decoupling of the impurity 
band from the main band, Vdec (in the Lifshitz model), or 
complete filling of all the states up to the impurity level Simp 
(in the Anderson model). 


tion parameter U, a crossing of these two can take place. 
By contrast, Va (c) in the Anderson model is practically 
constant: VA (0) ~ 2eimp- 

It is readily seen from Eq. U^that a considerable reduc¬ 
tion of the MIT bias can be reached by driving the impu¬ 
rity concentration c close enough to Ccr- This is an essen¬ 
tial advantage of the Anderson model scenario compared 
to that of the Lifshitz model, also taking into account 
that growing c simultaneously improves the thermal sta¬ 
bility of tuned MIT. The other practical advantage here 
is in a much higher steepness sm = d{£c — ep) /dV of 
this transition. This is seen from the comparison of cor¬ 
responding values: Vm ~ 0.34A, sm ~ 0.24 for the Lif¬ 
shitz model in Fig. [fl and Vm ~ 0.08A, sm ~ 0.42 for 
the Anderson model in Fig. [SI The latter advantage is 
even more enforced by the fact that the damping of Fermi 
quasiparticles (defining the Drude resistivity of metal) at 
V —>■ Vm varies slowly in the Lifshitz model, as in Eq. 
1241 but grows exponentially in the Anderson model, as 
in Eq. 1421 enabling here an extremely strong variation of 
the doped system resistivity near the tuned MIT. 

Summarizing, the doped and biased silicene presents a 
suitable opportunity for realization of practical electronic 
devices with tunable electric resistivity over a very broad 
scale, from normal metallic to fully insulating (possibly 
accompanied by respective optical, thermal, etc. effects), 
under rselatively weak bias. This regime can be opti¬ 
mized by a proper choice of impurity atoms, their loca¬ 
tion within the crystalline structure, and their concentra¬ 
tion. The comparative analysis of two common models 
for impurity perturbation on the host electronic spec¬ 
trum indicates the Anderson hybrid model (adequate for 
transition or rare earth impurities in silicene) to be more 
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Appendix 

Quasimomentum sums over the Brillouin zone (BZ) 
commonly result in certain analytic functions of other rel¬ 
evant arguments (such as position vectors, energy, etc.), 
and their calculation is done by passing from sum to in¬ 
tegral 

4 = — [ fkdk, (48) 

N ^ VBZ Jbz 

where vbz is the BZ volume. In the present case, the 2D 
BZ consists of two equilateral triangles with side y/SK, 
each centered in a nodal point, and the above integration 
of a function /k that decays fast enough away from the 
nodal points can be approximately done in the radial and 
angular variables ip over the circle of radius A: 

— [ /kdk Ri f ^d^ [ dipf {i, ip). (49) 

VBZ Jbz Jo Jo 

Among all the functions that have no such decay and 
so do not admit such approximation, we distinguish the 
important case of plane wave, /k = e**^ “, providing an 
exact result: 

= (50) 

k 

where the discrete Dirac delta is 6 (n) = 0 for any n 
joining two lattice sites (or two interstices) and 5 (0) = 1. 

Now we apply these techniques to the calculation of 
the basic interaction functions A„ and Eh in the approx¬ 
imation of non-perturbed spectrum. Starting from the 
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definition in the Lifshitz model, Eq. 1101 we obtain (at 
£2 < 


Now let us consider the standard deviation of the ran¬ 
dom scalar self-energy Epi 


A, 


Tj (e ± E/2) ^ e**"'" 

N ^£ 2 _]/ 2 / 4_^2 

A2 «“S2_y2/4_^2’ 


Tj{e±Vl2) 

A 2 


Kq (n/rs), 


the result in Eq. [TH In particular, its asymptotics at 
Tg defines the broadening of localized impurity level 
at c <C Co, by the criterion elm Si ^ 1 , where: 


1 roc 

ImSi = Im ^ , _ .2 - — “ Ai.y)rdr 

n^O ^ “ Ja 

Tirl cqA^ tt coA^ 

~ —1 -T ~ — In. — - --. 

a2 E (e - eioc) co V [e - eioc) 


Therefore the above criterion takes place at je — eioc\ ^ 
(cqA^/E) as indicated after Eq. 1301 

For the Anderson model, with the definition by Eq. |33 
we have (at > E^/4): 


F„ 


+ 


iUnr2 £+ ^COS ip 
« £2_y2/4_^2 

|_^e*k"(e-p^cos(/ 2 ) 

k 

2(e2-EV4) Y-e**"'"(e + Ccos(p) 

N 2^ £2_y2/4_^2 • 1 ) 


= S2 - Ep^ (53) 

restricting the definition of Ep in Eq. [37] to its first 
term, of the lowest linear order in impurity concentra¬ 
tion. Then we present: 

Ep ~ (c ^imp') ^ ^ CnEnE_n, 
n^O 

Ep Ri ie-eimp)~'^ ^ c„c„-FnE_„Fn/E_„/, (54) 
n,n '/0 

where the random numbers Cn of impurity occupation at 
nth interstice take the values 1 with probability c and 0 
with probability 1 — c and the non-renormalized denomi¬ 
nators e — £imp correspond to the adopted precision to the 
lowest order in c. Using this in Eq. |33| and taking into 
account that = c, CnCn' = c^ at n ^ n' and = c, 
we obtain: 

cr^ = c (1 - c) (e - eimp)~^ ^ [FnF^nf ■ (55) 

The decisive point in the evaluation of this sum over 
interstices is that the products FnF-n do not include 
Fourier components with coinciding quasimomenta (see 
after Eq. dni). Using the inverse orthogonality relation 
to Eq. [301 A^“^ = (5(k), we present the relevant 

sum as: 

Y^iFnF-nf = -F^ + S, 

n /0 


Then, applying Eq. [331 to the first sum in the right hand 
side of Eq. [51] gives two terms proportional to 6 (n) and 
dS (n) /dux, that vanish for any n 7 ^ 0. The second sum 
is treated by means of Eq. [331 as follows: 


2 

N 


E' 


.ik-] 


j e + C cos (fi 

£2 _ y2/4_^2 


2 

ttA^ 
e Jq 


/' 




£ 2 _E 2 / 4-^2 


( ilL 

\hvF 


(Jl 


(—) 

\hvFj 


cos 9 


and, after extending the upper integration limit to infin¬ 
ity, the exact formulas can be used: 


r°° xJq {x) dx 

Jq 52 - a ;2 

f°° x^Ji (x) dx 
Jq 62 - a :2 


ITT 

y 

iiT 


hP (b), 


6 ijf ’ ( 6 ), 


(52) 


with the 2nd kind Hankel functions {x) ( 2 ^ and 6 = 
n/|re|. This leads finally to the result of Eq. [331 


where Eq 


(cqe) is prevailed by 


s = E 


/ki /k2 /ka /ki-|-k2 -ka 


ki,k 2 ,ka 


with 


/k = 2 {s^ 


EV4) 


e -I- ^ cos (/? 

£ 2 _E 2 / 4 -^ 2 - 


This triple sum in kj is dominated by its real part, 
mainly due to the short-wave contributions by ^ ^ 
\J— E2/4 R! yTioA, and it can be estimated by the 
triple integral in = 6 .uf%: 


^imr) ^" ^imv 


/' 

J Si, 


d^3 


(co^imp) , 1 


6 (4^ + Cl + el) 


In — (56) 

Ccr 


(the (^-oscillating terms in /k are not important for this 
result). This readily leads to the expressions in Eqs. 1331 
and 1331 
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Finally, the broadening of localized impurity level at 
V 12 > £imp is obtained in the same way as shown above 
for the case of the Lifshitz model at c <C cq, with only 


difference for the pre-exponential factor cqA, here result¬ 
ing from the dominant contribution by the second term 
in Eq. [43l 
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